Universality of the rainfall phenomenon 
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We show that the universal properties of the rainfall phenomenon are the scaling properties of 
the probability density function of inter-drop intervals during quiescent periods, time intervals of 
sparse precipitation, and the universal shape of the probability density function of drop diameters 
during non-quiescent periods, time intervals of active precipitation. Our results indicate that the 
continuous flux-like vision of rainfall based on quantities such as the rain duration, rain intensity and 
drought duration is ineffective in detecting the universality of the phenomenon. A comprehensive 
understanding of rainfall behavior must rest on the acknowledgment of its discrete drop-like nature. 



INTRODUCTION 



Rainfall is a discrete drop- like phenomenon that has of- 
ten been described as a continuous flux-like phenomenon. 
The most common instrument used to measure the rain, 
the pluviometer, collects the water volume fallen through 
a given area per unit of time. The use of pluviometers 
and the importance of knowing intensity and duration of 
the rainfall phenomenon has lead to a description based 
on flux-like quantities such as the rain duration, rain in- 
tensity and drought duration (2j, even when radar mea- 
surement are used to infer the precipitable volume of rain 
[l3|. A flux-like view of rainfall is also central to the 
random cascade formalism used to describe rainfall pat- 
tern both in time and space [e.g., 0, [l(| EH- From 
a drop-like perspective, a considerable amount of work 
[e.g., [l], [H, 0, [ll| has been dedicated to the study the 
properties of raindrop spectra: the number of drops per 
diameter millimeter interval per cubic meter of air. Dou- 
ble stochastic Poisson processes have been used to de- 
scribe the variability in the drop counts per unit interval 
[e.g., M, EH- Only recently @, H[ an extensive study has 
been done of the properties of the sequences of inter drop 
time intervals and drop diameters as measured on the 
ground by disdrometers. Due to its importance for many 
aspects of human life, the rainfall phenomenon has been 
widely investigated. However, few works [e.g: |T3| are ex- 
plicitly dedicated to discuss the universal properties of 
the rainfall phenomenon. No entries were found (using 
common literature search engines) for works containing 
in their titles both the words "universality" (or univer- 
sal) and "rainfall" (or rain). So, are there any properties 
that a rain shower in New York and one in Rome have 
in common? 

In this letter, using data from the Joss Waldvogel im- 
pact disdrometer RD-69 located at Chilbolton (UK), we 
provide evidences that the universal properties of rain- 
fall phenomenon lie in the properties of drop-like quan- 
tities such as the inter drop time interval and the drop 
diameter. We show that 1) the flux- like view of the rain- 
fall phenomenon is not adequate to capture its universal 
features. 2) The temporal variability of rain can be de- 



scribed in terms of quiescent periods, periods of sparse 
precipitation characterized by a small drop diameter av- 
erage, and non quiescent periods, periods of active pre- 
cipitation characterized by a large drop diameter aver- 
age diameter. Moreover, the average and variance of the 
sequence of drop diameters are not stationary. 3) The 
probability density function of the inter drop time in- 
tervals t has an universal feature: a power law regime 
in the region r>lmin and r<lh. Inter drop time inter- 
vals in this range belong to quiescent periods. Finally, 
an universal shape for the probability density function 
of drop diameters during non-quiescent periods emerges 
upon removal of the non stationarity of the sequence of 
drop diameters. 



II. FLUX-LIKE VIEW OF A DROP-LIKE 
PHENOMENON 



The flux-like view of rainfall is that of a ON-OFF pro- 
cess. The rainfall time series is an alternating sequence of 
consecutive time intervals of duration A>0, the integra- 
tion time of the instrument used to monitor the precipi- 
tation, with (ON) or without (OFF) detectable precipi- 
tation. The relevant quantities [e.g.,[l3| are the duration 
of ON and OFF periods (rain duration and drought du- 
ration) and the volume of rain fallen during consecutive 
OFF periods (rain intensity). An integration time A>0 
causes all drop time intervals of duration r<A to be lost 
(equivalent to be detected as a drought of null duration) , 
and all inter drop time intervals of duration t>A to be 
detected as drought of duration [t/A] — 1 or [r/A] ([.] 
indicates the integer part). Thus, in the limit r^>A 



P£(r) cx A^(r) 



(1) 



where P&(t) is the distribution of drought durations and 
4>(t) is probability density function of inter drop time 
intervals. Eq. |T]) shows that distributions of drought 
durations relative to different integration times A will 
all have the same features at large times (t^A) as con- 
firmed by panel (a) of Fig. [T] This property is lost for 
the distributions of rain durations P&(t), panel (b) of 
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Fig. [TJ and rain intensities P^(v), panel (c) of Fig. [TJ As 
the integration time A increases a larger amount of inter 
drop time intervals are lost and rain durations that were 
separated by a drought duration are now detected as a 
longer rain period. This effect produces larger rain du- 
rations and, together with the temporal ordering of the 
sequence of drop diameters, larger rain intensities. 
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FIG. 1: Panel (a): Log-log plot of the probability Pa(0>t) 
of having a drought duration larger than r during the inter- 
val of time from January 24th to May 11th 2004. Different 
symbols indicate different values on the integration time in- 
tervals: A=10 seconds (white squares), A=l minute (white 
circles), A=10 minutes (white triangles), A=l hour (white di- 
amonds). Panel (b): Log-log plot of the probability Pa(#> t ) 
of having a rain duration larger than r during the same inter- 
val of time, and for the same the integration time intervals of 
panel (a). Panel (c): The probability P}±(V>v) of having a 
rain duration larger than v during the same interval of time, 
and for the same the integration time intervals of panel (a). 

The properties of the sequences of inter drop time in- 
tervals and drop diameters do not straightforward trans- 
late into those of the distributions of rain durations 
(P t a (t)) and intensities (P A (w)) [e.g., El]. Thus, these 
quantities are not valuable proxies to investigate the uni- 
versal properties of the rainfall phenomenon. 



However, some properties of the temporal variability of 
the rainfall phenomenon can be inferred. 



A. Quiescence 

Quiescence is a way of describing the temporal vari- 
ability of the rainfall phenomenon based on the relation- 
ship between inter drop time intervals and drop diame- 
ters. We divide the rainfall time series in time intervals 
of length S. A couple of consecutive time intervals of 
detectable precipitation may be separated by £;>0 con- 
secutive droughts (empty intervals) or be adjacent. In 
this case we say that the couple is separated by k=0 
consecutive droughts. Finally, for each couple we eval- 
uate the average number of drops n avg , and the aver- 
age drop diameter d avg . Fig. [2] shows the relationship 
between the couple average number of drops, the cou- 
ple average drop diameter and the number of consec- 
utive droughts in between the couple for (5= 10s. We 
see the tendency for a number of consecutive droughts 
k>l (interdrop time intervals r>10s) to separate cou- 
ples with a small average number of drops (<5) and a 
small average drop diameter (<0.6mm). To quantify this 
tendency, we introduce the concept of quiescence of order 
((5, to, n). A couple of consecutive time intervals of length 
5 with detectable precipitation and separated by k con- 
secutive droughts is a quiescent couple of order (5, m, n), 
if: 

k > m or k < m but rt avg < n, (2) 
where n avg is the couple average number of drops. 




III. TEMPORAL VARIABILITY 

Our data are from a Joss Waldvogel impact disdrom- 
eter RD-69 with a time integration A=10s, and 127 dif- 
ferent diameter classes (from 0.2998 mm to 4.9952 mm). 
Thus, we can record neither all the inter drop time inter- 
vals less than 10s, nor the exact arrival ordering of drops. 



FIG. 2: The logarithm of the density p(k,n avg ) of the couples 
(k,n avg ) in the kn a v g plane (bottom x-axis and left y-axis). 
The logarithm of the density p(k,d avg ) of the couples (k,d avg ) 
in the kd a v g plane (top x-axis and right y-axis). The in- 
ner plots show the logarithm of the densities p(k,n a v g ) and 
p(k,d avg ) in the region close to respective origins. Both den- 
sities are obtained using the data relative to the time interval 
from January 24th to May 11th 2004. 
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Quiescent (non-quiescent) periods are rainfall periods 
occupied by consecutive quiescent (non-quiescent) cou- 
ples. During quiescent periods, large inter drop inter- 
vals (r>m8) are preceded and followed by few (n avg <n) 
drops of small diameter (Fig. [2]). Non-quiescent or ac- 
tive periods are characterized by small inter drop time 
intervals (r<mS) separating drops with a larger range 
of diameters (Fig. [2]). Non quiescent periods are respon- 
sible for the bulk of precipitation: for a quiescence of 
order (5=10s, m=l, n=5) >95% of the precipitated vol- 
ume of water belongs to non quiescent periods. Some 
care is necessary in choosing the duration 8 of the time 
intervals dividing the rainfall time series: a 8 too large 
(>10 minutes) will result in a mixing quiescent with non- 
quiescent periods, a 8 too small (<1 second) will re- 
sult in too many time integration intervals with just one 
drop. Two quiescence of order (8,m,n) and (S',m',n') 
are "equivalent" if m'=m(8/8') and n'=n(8/8'). 



B. Drop diameter variability 

[f| introduce the concept of averaged "instant" shape 
to characterized the variability of raindrop spectra and 
their dep arture from the exponential form observed by 
fill ]. [Tif show that the probability density function r}(d) 
of drop diameters may change according to the portion 
(e.g.: dissipative edge, cloud base) or the type of storm 
(e.g.: orographic, non orographic) observed by a dis- 
drometer. Here, we show the temporal variability of the 
sequence of drop diameter is characterized by a moving 
average and a moving variance. Fig. indicates that the 
sequence of drop diameters does not have a constant av- 
erage, as quiescent periods have a lower average diameter 
than non-quiescent periods. A closer examination of the 
sequence of drop diameters reveals that its average to- 
gether with its variance are not stationary. Support for 
this thesis comes also from the results of [8[ . They report 
for the autocorrelation function of the sequence of drop 
diameters with a slow decay (the autocorrelation function 
reaches zero at lag «1250) followed by a long negative tail 
(lag >1250). The auto correlation function (as confirmed 
by simulations not reported here for brevity) of sequences 
of drop diameters exponentially and normally distributed 
with changing intensity around a moving average has the 
same features of that of m. 



IV. UNIVERSALITY 

In Fig. [3j we plot the probability P^(9>t) of having a 
drought duration larger than r for several time intervals 
of continuous observations at Chilbolton over a period of 
almost 2 years. All curves show a power law regime in the 
region between r>lmin and T<lh. The extensive period 
of time covered by our data, together with observations in 
other location of Earth's surface , indicate that the 

power law regime in the region from 1 minute to 1 hour of 



Fig. [3] is an universal property of the probability density 
function ip{ T ) of inter drop time intervals (Eq. [1]). This 
power law regime is a characteristic of quiescent periods: 
all quiescent couples (Eq.[5|) of order (8= 10s, m=l, n=5) 
are separated by inter drop time intervals r>10s (Fig. [2]). 
Moreover, Fig. [3] suggests that the end of the power law 
regime at r«lh signals a time scale separation between 
two different dynamics: the inter storm dynamics where 
quiescent and non-quiescent periods alternate each other, 
and the dynamics regulating the occurrence of different 
storms (meteorological dynamics). Thus, the probabil- 
ity density function ip(r) of inter drop time intervals can 
be thought as the sum of three components: 1) V'nq(t), 
the probability density function of non quiescent periods 
(tg[0s, < 10s]). 2) ?/>q(t), the probability density func- 
tion of quiescent periods (r€ [> 0s, < lh] ) with a power 
law regime in the region between T>lmin and r<lh. 3) 
ipQM (t) , the probability density function describing the 
meteorological variability (r>lh) of the particular loca- 
tion where the measurements are done. The index Q in 
tJjqm{t) indicates that all inter drop time intervals r>lh 
belongs to quiescent periods (Eq. [5] and Fig. [21). 
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FIG. 3: Log-log plot of the probability Pa(6>t) of having a 
drought of duration 6 larger than r as a function of r for dif- 
ferent time intervals of continuous observations at Chilbolton. 
The period of continuous observation were: from April 1st to 
November 3rd 2003 (white squares) , from November 5th 2003 
to January 5th 2004 (black squares), from January 8th to 
January 20th 2004 (black circles), from January 24th to May 
11th 2004 (white circles), from May 14th to July 17th 2004 
(white upward triangles), from July 19 th to August 2nd 2004 
(white downward triangles), from August 4th to August 19 th 
2004 (white diamonds), from December 10th 2004 to Febru- 
ary 28th 2005 (white pentagons). The curves are shifted for 
clarity. The ticks on the y axis indicates different decades. 

In panel (a) of Fig. |3J we plot the probability den- 
sity functions riq{d) of drop diameters of quiescent pe- 
riods relative to different months of observations. The 
observed variability is due to the non stationary char- 
acter of the sequence of drop diameters (Sec. IIII B|) . In 
fact, if the non stationarity is removed an universal shape 
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for the probability density function emerges. We con- 
sider non-overlapping time intervals of duration T and 
remove the average in every time interval. Panel (b) 
of Fig. 0] shows that the probability density function 
VQ(dT,[i) of the zero-average drop diameters sequence has 
a much smaller variability than the probability density 
function r)Q(d) of the original sequence (Fig. 2] panel (a)). 
If together with moving average also the moving vari- 
ance is eliminated (e.g. rescaling to unity the variance 
in each time interval), the probability density functions 
VQid-T.^.a) of the zero-average unitary- variance drop di- 
ameter sequences relative to different months "collapse" 
into a single curve (Fig. [4] panels (b) and (c)). The shape 
of this curve is not appreciably altered either by the 
choice of time intervals of different duration T (ranging 
from 10 seconds to ~10 minutes) to remove the non sta- 
tionarity of the sequence of drop diameters of quiescent 
periods, either by the use of different quiescence orders 
(n of Eq. @ ranging from 5 to 20) and of their "equiv- 
alence" classes (m and 8 of Eq. © changing in such 
a way to preserve the factors m6 and nS). The prob- 
ability density function rjQ^T^^) of the zero-average 
unitary-variance drop diameter sequence has two asymp- 
totic exponential tails: one for the positive and one for 
the negative values of the rescaled zero-averaged diame- 
ters (Fig. 2] panels (b) and (c)). A least squares fit of the 
exponential tails produce the following values for the de- 
cay constants: A+=2 (2mm<dr :A1(T <6mm) and A_=4.56 
(— 4mm<dj\ A1 . (T <— 2mm). 



V. CONCLUSIONS 

We introduce the concept of quiescence to describe the 
temporal variability of the rainfall phenomenon. The qui- 
escence captures a fundamental relationship (Fig. [2]) be- 
tween inter drop time intervals, drop diameter and their 
time ordering. These properties are not detected by the 
flux-like quantities such as rain duration and rain inten- 
sity. Using the concept of quiescence, we identify what 
are the universal properties of the rainfall phenomenon. 
The scaling property of the probability density func- 
tion of inter drop time intervals during quiescent periods 
(Fig. [3]) and the universal shape for probability density 
function of drop diameters (Fig. 2]). Our results sug- 
gest that the analysis of inter drop time intervals and 
drop diameters sequences and their properties offers a 
deeper insight than the analysis of the properties of flux- 
like quantities such as rain duration and rain intensity. A 
comprehensive understanding of the rainfall phenomenon 



must rest on its drop-like nature. 
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